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More Accurate Predictions with Transonic Navier-Stokes 
Methods Through Improved Turbulence Modeling 

by 

Dennis A. Johnson 

Experimental Fluid Dynamics Branch 
NASA Ames Research Center 


Because the aerodynamic characteristics of aircraft in the transonic 
regime are so sensitive to viscous effects, the selection of the turbulence model 
for a transonic prediction method is no less important than the selection of the 
numerical algorithm. Yet, the usual practice in transonic airfoil, Reynolds- 
Averaged, Navier-Stokes codes has been to employ "equilibrium" algebraic tur- 
bulence models. Satisfactory results are obtained with these turbulence models 
for weak interaction cases (i.e., cases where the upper surface shock wave is too 
weak to have a major effect on the turbulent boundary layer). Such is not the 
situation for cases where the shock wave is sufficiently strong to cause 
separation. The danger in using these "equilibrium" turbulence models for 
airfoil design is that they can result in unduely optimistic projections of aircraft 
performance at off-design conditions. 

Significant improvements in predictive accuracies for off-design con- 
ditions are achievable through better turbulence modeling; and, without neces- 
sarily adding any significant complication to the numerics. One well established 
fact about turbulence is it is slow to respond to changes in the mean strain Field. 
With the "equilibrium" algebraic turbulence models no attempt is made to model 
this characteristic and as a consequence these turbulence models exaggerate the 
turbulent boundary layer's ability to produce turbulent Reynolds shear stresses 
in regions of adverse pressure gradient. As a consequence, too little momentum 
loss within the boundary layer is predicted in the region of the shock wave and 
along the aft part of the airfoil where the surface pressure undergoes further 
increases. 

Recently, a "nonequilibrium" algebraic turbulence model was formulated 
which attempts to capture this important characteristic of turbulence. This 
"nonequilibrium" algebraic model employs an ordinary differential equation to 
model the slow response of the turbulence to changes in local flow conditions. In 
its originial form, there was some question as to whether this "nonequilibrium" 
model performed as well as the "equilibrium" models for weak interaction cases. 
However, this turbulence model has since been further improved wherein it now 
appears that this turbulence model performs at least as well as the "equilibrium" 
models for weak interaction cases and for strong interaction cases represents a 
very significant improvement. 

The performance of this turbulence model relative to popular 
"equilibrium" models is illustrated for three airfoil test cases of the 1987 AIAA 
Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 
"nonequilibrium" turbulence model is currently being applied to wing flows for 
which similar improvements in predictive accuray are being realized. 
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REYNOLDS SHEAR STRESS PROFILES 

AXISYMMETRIC BUMP. M m = 0.875 



197 





198 


X/C 




llT 

■j 

X § * 
0^0 
<2o 

F CC X 

So« 

— ^u. 

og° 

*gi 

ll t LU 
j □ DC 
< 2 h 

i o g> 

o z 

H Z O 
U. “J “ 

Ogo 

5 Q 

<§UJ 

: LU O 
CC SC 


— CO ^ 

LL (/) ? 

- Ir K 

H SF m 
< 0 - 5 : 
^ LU O 
0C CO ^ 
HI DC H 

fz > D 
Gqo 

CD < IE 


LU CD 

><o 

ffl £G 

Q DC 


199 


IMPROVEMENTS IN JOHNSON - KING TURBULENCE MODEL FOR 
WEAK SHOCK WAVE/BOUNDARY LAYER INTERACTIONS 

RAE 2822, M exp = 0.729, cv exp = 3.19° 
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